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Abstract 

The inelastic lifetimes of electrons in image-potential states at Cu(100) that are Stark-shifted by 
the electrostatic tip-sample interaction in the scanning tunneling microscope are calculated using 
the many-body GW approximation. The results demonstrate that in typical tunneling conditions 
the image state lifetimes are significantly reduced from their field- free values. The Stark-shift to 
higher energies increases the number of inelastic scattering channels that are available for decay, 
with field-induced changes in the image state wave function increasing the efficiency of the inelastic 
scattering through greater overlap with final state wave functions. 

PACS numbers: 73.20. At,68.37.Ef, 72. 15.Lh 
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INTRODUCTION 



The scanning tunneling microscope (STM) is a versatile and powerful probe of surface 
electronic structure; but it is not ideal. The electric field between the probe tip and the 
surface of the sample affects the surface. This influence can be exploited to positive effect, 
most dramatically through the controlled modification of surface atomic structure 
More prosaically, the influence of the tip must be allowed for when interpreting STM mea- 
surements, especially at semiconductor surfaces where tip-induced band bending occurs |j| . 

Recently a significant Stark-effect - the shift in energy due to the electric field - has been 
identified in scanning tunneling spectroscopy (STS) of Shockley surface-state electrons at a 
metal surface [4j. Surface states, in which the electron is caught between the barrier potential 
outside the surface and a band gap in the crystal, have been extensively studied by STM 
and STS, with a particular focus on their dynamics 

Baas 

and interactions [9], and the 

recognition of a Stark-shift in the case of Shockley states reconciles a discrepancy that has 
existed between STS-derived binding energies and those from photoelectron spectroscopy 
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FIG. 1: (a) Schematic energy diagram for the potential at the Cu(100) surface. The vacuum level 
lies near the center of the projected band gap which at ku = lies between the energies of the 
bulk X±i and X\ levels. The potential due to the image interaction (thick solid line) produces a 
hydrogenic-like series of excited states converging on the vacuum energy (dashed line); the n = 1 
level and corresponding probability density are shown, (b) Energy diagram for the surface in the 
presence of the tip of the scanning tunneling microscope, at the threshold bias for tunneling into 
the n = 1 image state. The field in the tunnel junction causes a Stark shift of the image state 
spectrum, with concomitant modification of the probability densities of the states. 
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12 1 in the case of a 



[lol ] . A more pronounced Stark-effect has been known for some time 
second class of surface electron state, namely the image-potential states that arise when an 
electron outside a conductor polarises the surface and is attracted to the resulting "image 
charge" , shown schematically in Fig. ^ Image-potential states are more weakly bound than 
Shockley surface states (which lie close to the Fermi energy Ep), forming a hydrogenic-like 
series with energies 

0.85 eV 

E n = E v - n=l,2,... 1 

(n + a) 2 

converging on the vacuum level of the surface E Y . In (0) a is a quantum defect that depends 
upon the surface. Tunneling via image states requires significantly greater bias voltages 



than Shockley states, and t 



several tenths of an eV 



re image state electrons are Stark-shifted to higher energies by 



The presence of a measurable Stark-shift in the surface state energies raises the important 
question as to whether there are also changes in the inelastic interactions of the surface state 
electrons in the presence of the STM tip. Electronic excitations in the surface state bands 
decay on a femtosecond timescale through interactions with the electrons and phonons of 
the surface and bulk, and there has been considerable activity in recent years directed at an 
understanding of these interactions A significant electric-field induced lifetime change 
would have important consequences for the interpretation of STM and STS experiments in- 
vestigating the dynamical properties of image-potential states, for example in nanostructures 
where the lateral resolution of the the STM is paramount. 

To investigate this issue we have performed many-body calculations of the lifetimes of 
image-potential states at Cu(100) in the presence of an electric-field due to the tip of an 
STM. Our calculations are based upon the approach introduced by Chulkov et al. [l5j and 
used subsequently in numerous surface state lifetime studies with considerable success 
The damping rate or inverse lifetime of the image state is calculated from the expectation 
value of the imaginary part of the non-local self energy operator 

r = h/r = -2 Jdr J dr 'ij? '(r)Im£(r , r'; E)ip(r'). (2) 

The energy of the image state is E, and ip the wave function. The imaginary part of the 
self energy is calculated in the GW approximation of many-body theory, which uses the first 
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term only in the series expansion of £ in terms of the screened Coulomb interaction W: 

1 f e 

Im£(r,r';e) = / de'lmG{r,r';e')lmW(r,r';e-e'). (3) 

^ Je f 

We use the zero'th order approximation to the Green function; in the spectral representation 

i 

where the ipi(r) are one-electron eigenfunctions with eigenenergies E± and 5 a positive in- 
finitesimal. The screened interaction is evaluated in the random phase approximation (RPA) 

W(r,r';u) = V(r-r') + J dr x j dr 2 V(r- ri ) 

xx°(r 1 ,r 2 - } u)W(r 2 ,r';u), (5) 

where V(r) is the bare Coulomb interaction and x°(r, r'; cu) is the density-density response 
function of the non-interacting electron system: 

2 r^ F 

X°(r,r';u>) = / de ImG(r, r'; e) 

J 

x[G(r : r';e+uj) + G*(r : r';e-uj)]. (6) 

This GW-RPA approach has been shown to give decay rates for image states at Cu surfaces 
ihat are within 1 meV of those found using the more complete GWr-TDLDA approximation 
16^ , in which exchange-correlation effects that are omitted in GW-RPA are included in both 
(J5J) and ©. Note that the phonon contribution to the decay rate of image-potential states 
at Cu(100) is < 0.5 meV and so safely ignored here |l7|. 

We calculate the Green function using a one-dimensional pseudopotential which by con- 
struction reproduces the Cu(100) bulk band edge energies and Xi, and the energies of 
the unoccupied surface resonance and first image state at the field-free surface, and which 
also accurately predicts the energies of the higher image states ^ | . In the direction parallel 
to the surface we assume parabolic dispersion with effective masses (m*) fitted to ab-initio 
band structures. To model the influence of the STM we follow Limot et al. 4| and include 
a linear potential due to the bias voltage between the STM tip and the sample, and modify 
the image potential to include the multiple images present in the tunnel junction geometry. 
Using this potential we are able to reproduce the sequence of Stark-shifted image state en- 
ergies and increments in tip-sample distance observed at Cu(100) in z(V) spectroscopy by 
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FIG. 2: Calculated decay rates (T) and lifetimes (r) of the Stark-shifted n = 1 image state at 
Cu(100) (filled circles). E\ is the image state energy in the presence of the electric field, and E v 
the field-free vacuum level of Cu(100). The significance of the dashed line is explained in the text. 



Wahl et al, who were also able to describe them using a model that omitted the effect 
of multiple-images. 

In Fig. |21 we show the calculated damping rates (inverse lifetimes) of the Stark-shifted 
n — 1 image potential state at Cu(100). Calculation parameters have been systematically 
varied to ensure that decay rates are converged to within 1 meV. At the field-free Cu(100) 
surface the image state lies at an energy of —0.57 eV relative to the vacuum energy, and we 
find T = 18±1 meV for this state, corresponding to a lifetime of r = 37±2 fs. This compares 
well with the lifetime r = 40 ± 6 fs measured using time-resolved two-photon photoemission 



(2PPE) 



19, 



and r = 38 fs found in previous GW calculations [16]. In the presence of 
the electric field due to the STM tip, the image state electrons are Stark-shifted to higher 
energies. At currents of 0.1-1 nA the n = 1 level is observed in z(V) spectroscopy at bias 
voltages near 4.7 V, corresponding to an energy E\ ~ E v + 0.1 eV [13]. For the results in 
Fig. 121 the tip-sample separation has been varied so that the bias voltage coincides with the 
n = 1 image state energy in the presence of the electric field. This situation corresponds 
to the onset of tunneling into the n = 1 level (Fig. E^b)). In these conditions the resulting 
image state energy increases linearly with the applied electric field, and we find that the 
decay rates also increase linearly. The rate of change is dT/dEi = 0.037 ± 0.002, so that 
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FIG. 3: Decay rates of electrons in the n = 1 image state band at Cu(100) as a function of lateral 
energy E(k\\)—Ei. Filled symbols correspond to Stark-shifted electrons with energy E\ = 4.7 eV at 
k\\ = 0; circles are calculated decay rates and triangles plus associated error bars are experimental 
decay rates from Ref. [l3j multipled by 2 to correct for an error in the phase relaxation length used 
in that work Open symbols are field-free results; circles are calculated decay rates, triangles 



plus associated error bars are decay rates from 2PPE experiments 



when the image state is Stark-shifted to E% = E v + 0.1 eV the lifetime is only 15 fs, a 
reduction of 60% from the field-free value. 

Recently the phase relaxation time of electrons in the n = 1 image state at Cu(100) 
has been studied using the STM by Wahl et al. [3], who measured the spatial decay of 
quantum interference patterns near steps. This technique measures the lifetime of electrons 
with non- vanishing momentum Hkn parallel to the surface, corresponding to energies above 
the image state band minimum: E(k\\) = E\ + h 2 k^/2m*. In Fig. El we compare calculated 
decay rates as a function of lateral energy with those reported in Ref. 13(. There it was 
concluded that the STM tip did not substantially alter the dynamical properties of the 
image-potential states, but subsequently an error has been recognised in the identification 
of the phase relaxation length that was used [21] so that the values displayed in Fig. El 
have been multiplied by 2 to correct for this. The calculations are performed for fields 
which give E\ = 4.7 eV. Also shown are calculated results for the field-free case, along with 
values from 2PPE measurements [2^ which correspond to this case. The overall agreement 
between theoretical and experimental lifetimes shown in Fig. El is very good, and confirms 
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FIG. 4: Probability density \ifi\ 2 of the n = 1 image state normal to the Cu(100) surface (z; 
Z{ m is the position of the image plane) and the negative imaginary part of the self energy — ImS 
evaluated for kn = and at the image state energy, shown for z = z' . Solid lines are for the 
unperturbed image state, E\ = E v — 0.57 eV, and dashed lines for the surface where the image 
state is Stark-shifted to an energy E\ = E v + 0.08 eV. 

the existance of a significant field-induced change in the inelastic lifetimes of the image state 
electrons. 

We now consider the origin of this effect. In Fig. 0]we compare the probability density 
\ip\ 2 and the imaginary part of the self-energy ImS at the Cu(100) surface calculated both 
with and without the electric field caused by the tip of the STM. The changes in \ip\ 2 show 
that accompanying the Stark-shift to higher energies is a redistribution of the weight of the 
surface state towards the metal surface, which increases the spatial overlap with the non- 
local self-energy. By displacing the image state electron towards the surface, the inelastic 
scattering channels are rendered more efficient. Calculating V using the wave function of 
the Stark-shifted state but the self energy of the field-free surface at the unperturbed image 
state energy accounts for approximately three-quarters of the full increase in T, as shown 
by the dashed line in Fig. 121 

The remaining change in the decay rate originates in the increase in the magnitude of 
ImS that can also be seen in Fig. 0] We find that calculating the self-energy using in 
(jnj) either the screened interaction W(r, r f ; e) of the field-free surface or of the surface in 
the presence of the electric field gives comparable results, i.e., the changes in the electron 
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wave functions caused by the electric field of the STM tip do not have a significant effect 
on the screening of the Coulomb interactions that dominate the inelastic scattering of the 
image state. Instead, the change in ImS is due to the increase in the number of final states 
into which the image state can decay, i.e. the number of states between Ep and the image 
state energy, which is Stark-shifted to higher energies by the tip-surface interaction of the 
STM. Thus the decreased lifetime of the Stark-shifted image-potential states results from an 
increase in the number of final states available for inelastic scattering along with increased 
efficiency of inelastic channels due to the greater spatial overlap of initial and final state 
wave functions. 

Given the magnitude of the tip-induced change in the lifetimes of image-potential states 
it is worthwhile to consider whether similar changes affect STM-derived lifetimes of Shockley 

surface states 011 El 3. 

Our calculations for Cu(lll) indicate that the effect is minor. 
In this case at feii = the Shockley state lies at Ep — 0.435 eV which means that under 
typical tunneling conditions the electrostatic tip-surface interaction gives rise to fields that 



ing into image-potential states, and 
. This causes only a minor change in 



are 5-10 times smaller than those present when tunne 
the resulting Stark-shift is much smaller: 10-15 meV (22], 
the number of final states that are available for decay. Furthermore, unlike image potential 
states which lie predominantly outside the surface, much of the Shockley state lies inside 
the metal, and is screened from the tip-induced electric field; there is a negligible change 
in the wave function penetration in the presence of the fields. Overall the electron-electron 
scattering decay rate which contributes two thirds of the total decay rate of ~ 21 meV 
changes by less than 5% in the electric field. Although it is safe to dismiss the effect for this 
particular case, it is clear that the tip-induced field will have an increasingly important effect 
on the lifetimes of Shockley states that lie further from Ep, especially at positive energies. 

To conclude, we have used the many-body GW-RPA method to calculate the inelastic 
lifetimes of electrons in image states at Cu(100) including the electric field due to the tip- 
surface interaction in the STM. We find that under typical tunneling conditions the lifetime 
of electrons in the n = 1 image state band is reduced by some 60% compared to in the absence 
of the STM tip. The Stark-shift to higher energies increases the number of inelastic decay 
channels that are available, whilst the the electric field moves the image state electrons closer 
to the metal surface, which significantly increases the efficiency of the scattering channels 
due to increased spatial overlap with final state wave functions. This tip-induced change in 
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electron lifetimes must be taken into account when using the STM to study the dynamical 
properties of higher-lying surface electron states. 
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